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ABSTRACT 



Context. We present results from a kinematical study of the gas in the nucleus of the active SO galaxy NGC 3998 obtained from 
archival HST/STIS long-slit spectra. 

Aims. We analyzed the emission lines profiles and derived the map of the gas velocity field. The observed velocity curves are consistent 
with gas in regular rotation around the galaxy's center. 

Methods. By modeling the surface brightness distribution and rotation curve of the Ha emission line we found that the observed 
kinematics of the circumnuclear gas can be accurately reproduced by adding to the stellar mass component a compact dark mass 
(black hole) of Mbh = 2.7^2q x 10^ Mq (uncertainties at a 2 cr level); the radius of its sphere of influence (^sph ~ 0".' 16) is well resolved 
at the HST resolution. 

Results. The BH mass estimate in NGC 3998 is in good agreement with both the Mbh - A^bul (with an upward scatter by a factor of 
~ 2) and Mbh - o" correlations (with a downward scatter by a factor of ~ 3 - 7, depending on the form adopted for the dependence of 
Mbh on cr). 

Conclusions. Although NGC 3998 cannot be considered as an outlier, its location with respect to the Mbh - cr relation conforms 
with the trend suggesting the presence of a connection between the residuals from the Mbh - cr correlation and the galaxy's effective 
radius. In fact, NGC 3998 has one of the smallest values of among the galaxies with measured Mbh (0.85 kpc) and it shows a 
negative residual. This suggests that a combination of both cr and is necessary to drive the correlations between Mbh and other 
bulge properties, an indication for the presence of a black holes "fundamental plane". 

Key words, black hole physics - galaxies: active - galaxies: bulges - galaxies: nuclei - galaxies: kinematics and dynamics 



1. Introduction 

It is generally believed that supermassive black holes (SMBHs; 
~ 10^ - IO'^Mq) are a common, if not universal, feature in 
the nucl ei of n earby galaxies. Since the discovery of quasars 
( Schmid lil963h . it has been suggested that active galactic nuclei 
(AGNs) are powered by mass accretion onto a SMBH (Salpeter 
[l964; Lvnden-Bell 1969). This belief, combined with the ob- 
served evolution of the space density of AGNs and the high 
incidence of low-luminosity AGN- like activity in the nucleus 
of nsM-by galaxies (Heckm anll 19801 R daoz et a l1ll995LlHo et al 



1997a b: Bra atz et alil 997: Barth et al. 1998, 1999: Nag ar et al 



2002) implies that a significant fraction of galaxies in the Local 




[M; 



Supports to these beliefs came from studies of the centers 
of nearby ea rly-type galaxies, which revealed that most con- 
tain SMBHs (iKormendv & GebhardJ200lLlMerritt & Fen-are's^ 
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1200 Ih and that the BH mass Mbh correlates with some 
properties of t he host galaxy, such as bulge lumi nosity Lh,,! 
and mass Mh ,.i (iKo rmendy & Richstone .1995: Magorrian et alJ 
1998; Marcon i & Huni2003h . light concent ration (iGraham et all 
2001) and bulge velocity dispersion o-f,,,] ("Ferrarese & Merritj 
2000; Gebhardt et al. 2000; Tremaine et al. 2002). The coiTela- 
tion with the bulge velocity dispersion was thought to be the 
tightest having the smalle st scatter: rms ~ 0.3 in log Mbh- 
Ferrarese & Merrit^ ( l2000h argued that, for their selected sam- 
ple of 12 galaxies, thought to have the most reliable BH mass 
estimates, the observed scatter in the Mbh - crbui relation was 
fully accounted for by the assumed measurement errors, which 
implies that there may be no intrinsic scatter in the correlation. 
More recently, Marconi & Hunt (2003) have shown that when 
considering only galaxies with secure BH mass and bulge pa- 
rameters determinations a ll the above correlations have a similar 
observed scatter (see alsolMcLure & Dunlotl2002t lErwin et alJ 
l2()04llHarin'g & Rixl2004) . 

The existence of any intrinsic correlations of Mbh and 
host galaxy bulge properties supports the idea that the growth 
of SMBHs and the formation of bulges are closely linked 
(Silk&Rees 1998; Haehneh & Kauff^mann 2000), therefore 
having important implications for theories of galaxy formation 
and evolution. Moreover, SMBH mass estimates inferred via the 
above correlations, when more direct methods are unplayable, 
enter in a variety of important studies spanning from AGNs 
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physics and fueling to coeval formation and evolution of the host 
galaxy and its nuclear black hole. These correlations need, how- 
ever, to be further investigated by increasing the number of ac- 
curate BH mass determinations in nearby galactic nuclei to set 
these correlations on stronger statistical basis. In particular such 
a study has the potential to establish the precise role of the vari- 
ous host galaxy's parameters in setting the resulting BH mass. To 
date, reliable SMBH detections have been obtained for a limited 
number of galaxies (~ 30, see^errarese & Ford 2005 for a re- 
view), with a bulk of Mbh estimates in the range of 10^ - IO^Mq. 
Furthermore, the coverage of the mass range is particularly poor 
at the low-mass end, the sample being strongly biased against 
late Hubble types and low luminosity objects. To add reliable 
new points to the Mbh- host galaxy properties planes is then a 
fundamental task for future developments of astronomical and 
physical studies, and is the aim of the work described in this 
paper 

Spectral information at the highest possible angular resolu- 
tion is required to directly measure the mass of SMBHs: the 
radius of the"sphere of influence" ( iBahcall & Wolj ri976l) of 
massive BHs is typically ~ 1" even in the nearest galax- 
ies. Among the so far most widely used techniques to de- 
tect and estimate BHs masses is stellar dynamical modeling 
(e.g. Dressier & Richstone 1988; Kormendv & Richston e 1995.: 
van der Marel et al.lll99'8r iGebhardt et al. 2000; Verol me et all 
2002 1). but the interpretation of the data is complex involving 
many degrees of freedom an d requiring data of very high signal- 
to-noise ratio (IValluri et alj '2004). Radio frequency measure- 
ments of H2O masers in disks around BHs, finally, can be ap- 
pUed only to the small fraction of the disk s inclined such tha t 
their maser emission is directed toward us (iBraatz et al."1997'). 
A more widely applicable and relatively simple method to de- 
tect BHs is based on gas kinematics (e.g . Harms et al. 1994; 
iFerrarese et al]ll996t iMacchetto et"al]ll997: .Barth et aL.2001i) . 
through studies of ordinary optical emission lines from circum- 
nuclear gas disks, provided that the gas velocity field is not dom- 
inated by non gravitational motions. A successful modeling of 
the gas velocity field under the sole influence of the stellar and 
black hole potential is needed to provide a posteriori support 
for a purely gravitational kinematics. The angular resolution of 
~ (X' 1 of the Space Telescope Imaging Spectrograph (STIS) on- 
board HST is the most suitable to perform such studies. The 
wealth of unpublished data contained in STIS archives repre- 
sents an extraordinary and still unexplored resource. 

Following this chain of reasoning, we performed a system- 
atical search for unpublished data in the HST STIS archive with 
the aim of finding galaxies candidates to provide a successful 
SMBH mass measurement. 

In this paper we present the results obtained for the ac- 
tive galaxy NGC 3998. From the Lyon/Meudon Extragalactic 
Database (HyperLeda), NGC 3998 is classified as an early type 
(SO) galaxy with an heliocentric radial velocity of 1040 + 18 
km s"'. With Hq - 75 km s"' Mpc"' and after correction for 
Local Group infall onto Virgo, this corresponds to a distance of 
17 Mpc and a scale of 83 pc arcsec"' . NGC 3998 is spectroscop- 
ically classified as a LINER with broad Ha emission ( Heckman 
[1980; Keel 1983; Hoetal. 1997a b) and no significant broad- 
line polarization (Barth et aL.199 9) therefore indicating a proba- 
ble direct view to a BLR. Support to this hypothesis comes from 
studies of HST WFPC2 optical images which revealed an unob- 
scured nucleus and the presence of a bright circumnuclear ion- 
ized gas disk ( Pogge et al. 2000). The nucleus contains a var iable 
compact flat spectrum radio source jHummel et alJlT98l dis- 



playing a weak jet-like northern structure jFilho et ai]l2002h . A 
bright UV source (unresolved at HST FOC resolution) is present 
in the center of NGC 3998 (Fabbiano et al. 1994). Its rapidly 
variable flux ( Maoz et al. 2005 ) implies that a significant fraction 
of the UV output is contributed by a non-stellar AGN compo- 
nent. A substantial featu reless continuum coniponen t is also ob- 
served in the optical (Gonzalez Delgado et al. 2004). The pres- 
ence of an AGN in the nucleus of NGC 3998 is further strength- 
ened by X-ray observations, which reveal ed the presence of a 
nuclear source with a power law spectrum (iRoberts & Warwick 
2000; Pellegrini et al. 2000) and X-ray and Ha luminosity values 
consistent with the low luminosities extension of the correlation 
observed for Seyferts and QSOs (see K oratkar et al. 1995 ). 

The paper is organized as follows: in Sec. |2 we present 
HST/STIS data and the reduction that lead to the results de- 
scribed in Sec. 13 In Sec.|?]we model the observed Ha emission 
line rotation curve and we show that the dynamics of the circum- 
nuclear gas can be accurately reproduced by circular motions in 
a disk when a point-like dark mass is added to the stellar poten- 
tial. Our results are discussed in Sec.|5l and summarized in Sec. 

m 



2. HST data and reduction 

NGC 3998 was observed with STIS on HST in 1997 November 
01 with the G750M grating and the 52"x Of.' I sHt. Data were ac- 
quired at five different slit positions, following a perpendicular- 
to-slit pattern with a step of Of.' 1 and the central slit centered on 
the nucleus. The orientation of the slit was -81° from North, 
and the exposure time was 328 s for each position. The position 
of the continuum peak had been acquired with the ACQ mode 
with two 10 s exposures obtained with the optical long-pass fil- 
ter MIRVIS. The five spectra obtained, NUC for the nuclear slit, 
N1-N2 and S1-S2 (from North to South) for the four off-nuclear, 
were retrieved from the public archive. Fig.^shows the slit loca- 
tions superposed onto a narrow band WFPC2, obtained by using 
the F658N filter, including the Ha+[N II] lines. 

The data were obtained with a 2x1 on-chip binning of the 
detector pixels and automatically processed through the standard 
CALSTIS pipeline to perform the steps of bias and dark subtrac- 
tion, applying the flat field and combining the two sub-exposures 
to reject cosmic-ray events. The data were then wavelength and 
flux calibrated with conversion to heliocentric wavelengths and 
absolute flux units and rectified for the geometric distortions. 
The 2-D spectral image obtained for each slit position has a spa- 
tial scale of 0'.'0507 pixeL' along the slit, a dispersion of A/l = 
1.108 A pixel ' and a spectral resolution ofR^ 3000, covering 
the rest frame wavelength range 6480-7050 A. 

For each spectrum we selected the regions containing the 
lines of interest. The lines were fitted, row by row, along the 
dispersion direction, together with a linear continuum, with 
Gaussian functions using the task SPECFIT in STSDAS/IRAF 
All emission lines present in the spectra (Ha, [NII]/l/}6548,6583 
and [SII]/l/l6716,6731) were fitted simultaneously with the same 
velocity and width and with the relative flux of the [Nil] lines 
kept fixed to 0.334. Results of the fit are tabulated in Appendix 

El 

A single Gaussian function for Ha line does not produce an 
accurate fit to the lines profile in the region r < Of.'l for three po- 
sitions of the slit: NUC and the nearest Nl and S 1 . Here the lines 
profile shows a broad base component superposed to the narrow 
lines of Ha and [N II] (see Fig.|2]for NUC slit), indicative of the 
presence of a Broad Line Region. In these regions we performed 
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Fig. 2. Nuclear spectrum of NGC 3998 showing the presence of 
a broad Ha component. The best model spectrum (dashed line) 
is superposed to the data; dotted lines are for the gaussians and 
linear continuum components of the fit. The residuals from the 
fit are shown in the lower panel. 



a fit with two Gaussians components for Ha line, a narrow and 
a broad one whose velocity we kept fixed at all locations. 



In the external regions, where the SNR was insufficient the 
fitting was improved by co-adding two or more pixels along the 
slit direction. 



3. Results 

The results obtained from the fitting procedure at the five slit 
positions are shown in Fig.|3]through Fig.|5]where we show the 
central velocity, surface brightness and FWHM for the narrow 
Ha line at each location along the slits. Emission is detected out 
to a radius of ~ (X'S corresponding to ~ 40 pc. The extension of 
the line emission and its behaviour along the five slit positions 
are in agreement with the Ha image seen in Fig.Q where an ~ 
EST- WEST elongated emission structure is visible. Conversely, 
the stellar component shows an almost circular distribution in 
the V band continuum. 

The line emission is strongly peaked on the central slit (Fig. 
13 middle panel), rapidly decreasing at larger radii and not show- 
ing the presence of any emission line knot other than the central 
maximum. The velocity curve in the central slit, NUC, has a 
full amplitude of ~ 400 km s 'and shows a general reflection 
symmetry: starting from the center the velocity rapidly rises on 
both sides by ~ 200 km s 'reaching a peak at r ~ O'.'l from the 
center. At larger radii, the velocity decreases to form a plateau 
before rising again at the extremes of the velocity field. Both the 
line flux and the line width rapidly decrease from the nucleus 
outwards. 

The behaviour seen in the off'-nuclear slits is qualitatively 
similar to that seen at the NUC location, but with substantially 
smaller velocity amplitude and, more important, a less extreme 
nuclear gradient. Both amplitude and gradient decrease at in- 
creasing distance of the slit center from the nucleus, with a be- 
haviour characteristic of gas rotating in a circumnuclear disk. 
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Fig. 4. Same as Fig.|3for the nearest nucleus slit positions Nl (left panel) and SI (right). In the region r < 0'.'2 fits are performed 
with a double Gaussian for Ha line to separate the broad and narrow components. 
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Fig. 5. Same as Fig.|3lfor the N2 (left panel) and S2 (right) sht. 



4. Modeling the rotation curves 

It therefore appears that an ionized gas system is present in the 
innermost regions of NGC 3998 with a smooth and regular ve- 
locity fie ld, co-rotating with respect to the larger scale stellar and 
gas disk jFisheJl997i) . 



Ou r modeling code, described in detail in iMarconi et alJ 
()2003l), was used to fit the observed rotation curves. The code 
computes the rotation curves of the gas assuming that the gas 
is rotating in circular orbits within a thin disk in the galaxy po- 
tential. The gravitational potential has two components: the stel- 
lar potential (determined in the next section), characterized by 
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Fig. 3. Velocity, surface brightness and FWHM for the narrow 
Ha component along the NUC slit. Surface brightness is in units 
of 10"'^ erg s"' cm"^ arcsec"^. Positions along the slit are rela- 
tive to the continuum peak, positive values are NW. In the region 
r < (y.'2 fits are performed with a double Gaussian for Ha line to 
separate the broad and narrow components. 



Fig. 6. Results of the isophotes analysis of the V band image 
of NGC 3998. Surface brightness is shown in the top panel (in 



units of erg s" cm A" arcsec" ), the galaxy's ellipticity and 
position angle are shown in the middle and bottom panels re- 
spectively. 



its mass-to-light ratio and a dark mass concentration (the black 
hole), spatially unresolved at HST-t-STIS resolution and charac- 
terized by its total mass Mbh- In computing the rotation curves 
we take into account the finite spatial resolution of HST+STIS, 
the line surface brightness distribution and we integrate over the 
slit and pixel area. The is minimized to determi ne the free 
param eters using the downhill simplex algorithm bv lPress et alJ 

4. 1 . The stellar mass distribution 

In order to assess the contribution of stars to the gravitational 
potential in the nuclear region, we derived the stellar luminosity 
density from the observed surface brightness distribution. 

We reconstructed the galaxy light profile using a WFPC2 
F547M (V band) image retrieved from the public archive (Fig. 

An inspection of HST archive images shows that contamina- 
tion from dust is negligible, since dust absorption is seen only on 
relatively large scale (r > 2 kpc) w hile the cent ral regions are ap- 
parently free from dust (see also Pogge et al. 2000). We used the 
IRAF/STSDAS program ELLIPSE'to fit elliptical isophotes to 
the galaxy (see Fig.|6}. Excluding the nuclear regions (r < Ol'IS) 
that are dominated by a central compact source, the ellipticity 
shows small variations around a value of ~ 0.15. The position 
angle does not show significant variations being approximately 
constant at PA = 135°. 

The nature of the compact nuclear source is crucial for the 
estimate of the stellar mass distribution. In fact, if this is asso- 
ciated to light produced by the active nucleus, it does not corre- 
spond to a stellar mass contribution and should not be included 
in the mass budget. The presence of a bright (and variable) UV 



and X-ray source seems to indicate that indeed this is the case. 
We explored in more detail this issue by fitting the brightness 
profile with a Sersic law ( Sersic 1968T with superposed a point 
source, whose profile was derived from a synthetic Point Spread 
Function modeled with TINYTIM. This analysis shows a very 
good agreement between the data and the model (see Fig. 0. 
This supports the conclusion that the central source is unresolved 
and associated to AGN emission. 

The inversion procedure to derive the stars distribution from 
the surface brightness is not unique if the gravitational potential 
does not have a spherical symmetry. Assuming that the gravita- 
tional potential is an oblate spheroid, the inversion depends on 
the knowledge of the potential axial ratio q, and the inclination 
of its principal plane with respect to the line of sight. As these 
two quantities are related by the observed isophote ellipticity, 
we are left with the freedom of assuming different galaxy in- 
clinations to the line of sight. We performed the de-projection 
adopting for the galaxy inclination the value / - 40°, as given 
from the HyperLeda database. Due to the small ellipticity of this 
galaxy, the precise value of its inclination has only a marginal 
effect on the resulting mass distribution. 

Following van der M arel & van den BoschI Jl998h . we as- 
sumed an oblate spheroid density distribution parameterized as: 

where m is given by = + + z^/q^, xyz is a reference 
system with the xy plane corresponding to the principal plane 
of the potential and q is the intrinsic axial ratio, and performed 
the de-projection adopting for the mass-to-light ratio in the V 
band the reference value of Ty - 1. A detailed description of 
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Fig. 7. Fit to the brightness profile of NGC 3998 obtained with 
a Sersic law with an added nuclear point source (dashed line). 
The residuals from the fit are shown in the lower panel. 

the relevant formulas and of the inversion and fit procedure is 
presented in'Marconi et al.' ('2003V The best fit obtained is shown 
in Fig.|S]with a = 1.67, yS=0.65 and r\, =4.65"'. A point source 
with flux 5.2 10"'* erg s"' cm"- A"' was added to the extended 
luminosity distribution. 

4.2. Fitting tine gas kinematics 

Our modeling code was used to fit the nuclear rotation curves. 
The free parameters of the fit are: 

- the systemic velocity, v^ys, 

- the impact parameter (i.e. the distance between the nuclear 
slit center and the center of rotation) b, 

- the position of the galaxy center along the nuclear slit so, 

- the angle between the slits and the line of nodes, 0, 

- the disk inclination /, 

- the mass-to-light ratio of the stellar component, Yy, 

- the black hole mass Mbh- 

In a oblate spheroid, the stable orbits of the gas are copla- 
nar with the principal plane of the potential and it is possible to 
directly associate the galaxy inclination and line of nodes with 
those of the circumnuclear gas. However, the potential shape is 
not sufficiently well determined by the isophotal fitting down 
to the innermost regions of the galaxy and it is possible that a 
change of principal plane might occur at the smallest radii, in 
particular within the sphere of influence of a supermassive black 
hole. By these considerations we preferred to leave the disk in- 
clination as a free parameter of the fit. We then performed a 

' Note that with the derived density profile, the stellar mass included 
within the HST spatial resolution (UH) is 4.0 xIO'Mq, having adopted 
a mass-to-light ratio Ty = 6.5 that will be derived in the next Section. 
Thus it only represents a fraction of ~ 15 % of the best fit value for the 
SMBH mass. 



Fig. 8. Fit to the surface brightness profile obtained from an 
oblate spheroid stellar density distribution with an added nuclear 
point source. 

minimization for different values of /, namely / = 10°, 80°, al- 
lowing all other parameters to vary freely. Due to the sensitivity 
of the observed line width to the brightness distribution model- 
ing and to other co mputational problems (i.e. a coarse sampling, 
iMarconi et al 72006) we decided to initially perform the fit with- 
out using the values of the line width. We will show in Sect. 
I4.2.1l that the inclusion of the line widths in the modeling code 
has only a marginal effect on our results. 

To build the synthetic kinematic al models the intrinsic line 
surface brightness distribution for the narrow Ha line had to be 
obtained for each disk inclination. The observed emission line 
surface brightness was modeled for each / value with a composi- 
tion of three circularly symmetric Gaussian functions. The first 
reproducing the observed central emission peak, the second the 
intermediate regions and the third to account for the brightness 
behaviour at large radii. The modeling was performed through a 

minimization, leaving as free parameters the intensity, scale 
radius and peak position. The choice of a particular model for the 
line surface brightness distribution does not affect the final BH 
mass estimate provided that the model reproduces the observed 
line emission within the errors. Nevertheless it has an important 
effect on the quality of the velocity fit (Marconi et al. 2006). 

The x^ reduced values (x^ - x^/d.o.f.) of the best fit to the 
velocity curves, obtained for each disk inclination, are shown in 
Fig. |5] The quality of the fit depends only very weakly on the 
assumed disk inclination, with almost constant values of Xr up 
to / - 60°, while rapidly increases at larger inclinations. The 
overall best fitting model to our data is obtained for / - 30° for 
the set of parameters reported in Tableland is presented in Fig. 

M 

The value of minimum is far larger than the value indica- 
tive of a good fit and this is in contrast with the fact that the 
curves shown in Fig. ^| seem to trace the data points well. The 
reason for this discrepancy is that;t'j^ is not properly normalized 
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< 70° yield acceptable fits (within the 2cr level), apparently 
^ 1 our analysis only provides a lower bound to the black-hole mass. 

However, also the mass-to-light ratio Ty increases sharply at low 
inclinations reaching a value of 50 for / = 10°. In fact at smaller 
inclinations a deeper potential is needed to maintain the pro- 
jected velocities seen at the largest radii where the stellar compo- 
nent dominates. We then considered the evolutionary synthesis 
models for stellar populations derived by Maraston ( 1998) in or- 
9 h / - der to limit ourselves only to astrophysically acceptable models. 

The value of Ty monotonically increases with the galaxy's age 
reaching, for a Salpeter initial mass function (IMF), the maxi- 
mum value in the V band of Ty - 7.87 for an age of 15 Gyr 
This value corresponds to an inclination of ~ 27° (having used a 
1/sin^ i interpolation on the disk inclination for the dependence 
of both Ty and Mbh)- This translates into an upper bound to the 
I black hole mass (see Fig. II 11 1. Therefore, although larger black- 
hole masses are allowed when considering only the kinematical 
modeling, they are in reality unacceptable as they would corre- 
spond to a mass-to-Ught ratio of the stellar population larger than 
predicted by the stellar population synthesis. Furthermore, the 
J — J adopted upper limit on Ty compares favourably with observa- 
tions. The corr elation be tween galaxy' s luminosity an d mass-to- 
1 Ldegi eej yight ratio (e.g. lvan derMarel.i99L.Cappellari et alj.200 6) pre- 

dicts a range (in the R band) of Tj? ~ 3 - 6.5. This is consistent 
Fig. 9. xi values of the best fit obtained for each disk inclination, with an age of 7 - 15 Gyr and with the upper limit we adopted 

in the V band. 
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Table 1. Best fit parameter set 



/ b 


^0 


e 




Tv 


MBH(Afo) 


xi 


30 0.00 


-0.02 


28 


1010 


6.5 


2.7x10"* 


5.96 



By combining the allowed range in inclination with the sta- 
tistical uncertainties on Mbh associated to variations in Ty, we 
obtain a global range of acceptable black hole mass of Mbh — 
(0.7 - 5.1) X lO'^Mo at a 2 0- level. 



(e.g. because not all points are independent or as they do not in- 
clude the uncertainties in the relative wavelength calibration for 
the five slits) and/or imply t he presence of s mall deviations from 
pure rotation. Following Barth et alJ J2001h we then rescaled the 
error bars in our velocity measurements by adding in quadrature 
a constant error such that the overall best fitting model provides 
X^ld.o.f. ~ 1. This is a quite conservative approach as it has the 
effect of increasing the final uncertainty on Mbh- The additional 
velocity error is found to be 28 km s We rescaled all values 
of with this procedure (Fig. II II bottom panel). The best fit- 
ting models obtained at varying disk inclination are within the 
2cr confidence level {^xl < 0.09) for 10° < / < 70°. 

To evaluate the statistical uncertainty associated to the black 
hole mass estimate we explored its variation with respect to the 
parameters that are more strongly coupled to it, i.e. the mass- 
to-light ratio Yy and the gas disk inclination /. The uncertainty 
on Mbh associated to changes in Yy has been estimated build- 
ing ax^ grid in the Mbh vs. Yy parameter space. At each point 
of the grid, described by a fixed pair of Mbh and Yy values, 
we obtained the best fit model allowing all other parameters to 
vary freely and derived the corresponding Xr value (properly 
rescaled). This enabled us to build contours of confidence level. 
The result of this analysis at ; - 30° is presented in Fig.^] The 
2cr ranges of the BH mass and mass-to-light ratio at this disk 
inclination are Mbh = 2.7+[ q x 10*^ Mq and Yy = 6.5t_l l respec- 
tively. The statistical uncertainty on Mbh is reported in the error 
bar in Fig. We repeated the same analysis for another rep- 
resentative value of the inclination, 60°, obtaining similar frac- 
tional uncertainties on the parameters. 

The dependence of BH mass on the gas disk inclination 
closely follows the expected scaling with a oc 1/sin^ / law (see 
Fig. II 11 1, at least up to / < 60°. Since all disk inclination 



4.2.1. Line width distribution 

At this point of our analysis we tested the influence over the 
above results of including the line width in the fitting procedure. 
Adopting the best fit parameter set (Table ^ we obtained the 
FWHM distribution shown in Fig^] (solid line). The observed 
line widths are acceptably well reproduced by the model. Only a 
slightly underestimate (~ 20%) of the nuclear increase is observ- 
able at NUC, while the largest deviation from the peak values 
occurs at the off-nuclear S 1 . Among the line parameters, the line 
width is the most sensitive to the brightness distribution mod- 
eling and to other computational problems , i.e. a sub-samplin g 
of the grid used by the numerical code (Marconi et al. "200^. 
Due to this sensitivity only weak constraints on the BH mass es- 
timate can be derived by its inclusion in the fitting procedure. 
Therefore, a slightly underestimate of the peak observed values 
does not invalidate the derived fitting model. 

Despite the goodness of the result, we repeated the x^ min- 
imization at / - 30°, this time by including the observed line 
widths in the modeling procedure. The best fit obtained is shown 
in Fig. ^] (dashed line). Only a small improvement of the 
match with observed values is obtained, while the model velocity 
curves remain substantially unchanged respect to the previous 
result. Furthermore, BH mass and mass-to-light ratio values do 
not change significatively, resulting Mbh = 3.1 x 10^ Mq and Yy 
- 5.6. The modeling code assumption that the nuclear gas is in a 
thin, circularly rotating disk is then verified through the satisfac- 
tory good match of observed and model line width distributions. 
Indeed, the nuclear rise of the line width is well accounted for as 
unresolved rotation by the fitting model. 
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Fig. 10. Overall best fit to the rotation curves (left), assumed line surface brightness distribution (right). From upper to bottom 
panel: N2, Nl, NUC, SI and S2. 



5. Discussion 

Our model fitting of the nuclear rotation curves of NGC 3998 
indicates that the kinematics of gas in its innermost regions can 
be successfully accounted for by circular motions in a thin disk 
when a point-like dark mass (presumably a supermassive black 
hole) of Mbh - 2.7+2 Q ^ 10** Mq is added to the galaxy potential. 

Let us explore how this mass determination is connected 
with the properties of the host galaxy. We compare our BH mass 
estim ate with the k nown correlations with host spheroid (bulge) 
mass ( lMarconi"&^ Hunt 2003) an d with the stellar vel ocity dis- 
persion (Tremaine et al. 2002; Ferrarese & Ford 2005). 

Following Marconi & Hunt (200 3), we used the virial mass 
(Mvii = 3ReCrljG) to determine the bulge mass of NGC 3998. 

The effective radius of the bulge, R^, has been estimated 
as the weighted average of the determinations by Fisher et all 
( fT996l) (11") and[Sanchez-Portal et al. (2004) (9'.'6), from which 
we derived = 10'.'2 + ^.'6 (0.85 + 0.05 kpc). 

We found three optical determinations of the stellar ve- 
locity dispersion for NGC 3 998 i n the literature. The val- 
ues are 297 kms^'lFisheJj 119971 2"x 4" slit), 314 + 
20 kms"NTonrv& Davis 1981, 3"x 12" sht), 333 + 22 
km s''( Nelson & Whittleiil995i 1"5 x 2'.'2 sUt) respectively ^ 
We adopted the value from HyperLeda database o-jtar = (305 + 
10) km s"'. The velocity dispersions in the catalogue are mean 
values standardized to a circular aperture of radius r^p - 0.595 
h-' kpc (see lGolev & Prugnieli ri998). At the distance of NGC 
3998 this radius corresponds to an aperture of 9'.'6 of radius, very 
close to the value of the bulge effective radius for NGC 3998. 

Assuming that this value of iTstar is a good approximation 
for cTbui, we obtained for the bulge mass the value Mbui - 
(5.5 + 0.7) X 10'"Mq. Using this estimate of Mbui and the corre- 

2 lBertolaetalJ<1984) measured cr = 350 ± 40 km s ' , but, because 
of saturation, their data gave no resuUs on stellar kinematics in the in- 
nermost central regions (< 2"). 



lation o f^ arconi & HuntI (|200^ which considers only "secure" 
BH mass determinations (i.e. BH with resolved sphere of in- 
fluence: 2Rsph/Ries > 1, with /?res the spatial resolution of the 
observations) the expected Mbh for NGC 3998 is 1.3 xIO^^Mq, 
in excellent agreement, within a factor of 2, with our determina- 
tion (see Figll4>. We also note that the value of the BH sphere 
of influence radius in NGC 3998 (Rsph - GMbh/o-I^^i) is Rsph ~ 
13 pc (0'.'16). This implies a well resolved BH sphere of influ- 
ence at the HST resolution (~ O'.'l) for our Mbh determination 

(2^sph/^res - 3.2). 

Concerning the correlations between BH mass and central 
stellar velocity dispersion, ad opting the correlation parameters 
estimated by Tre maine et al. ' (2002) the expected BH mass for 
NGC 3998 is 7.4 xIO^Mq (see Fig[l5}, a factor ~ 2.7 higher 
than our estimate. However, NGC 3998 cannot be considered as 
an outlier from the Mbh vs. cr correlation given the errors in the 
estimates of Mbh, ostdi-, and of the best fit parameters describing 
Mbh - cr relation, and its intrinsic scatter (0.25-0.3 in log Mbh)- 
We can also compare our measurement with the expectations 
of the Mbh - o" relation recently derived by Ferrarese & Ford! 
(i2005V We then normalized the central velocity dispersion to an 
aperture of radi us equal to 1/8 of Re, following the method intro- 
duced by Jor gensen et alJ d 19951) . and derived ctr jg = 330 + 11 
kms This form of the correlation predicts a BH mass of 1.9 
xIO^Mq, a factor of 7 higher than our estimate, a larger discrep- 
ancy than the one found adopting the Tremaine et al. values. 

|m arconi (feHund (l2003h showed with a partial correlation 
analysis that Mbh is separately significantly correlated both with 
cr and Re- This is clearly shown by the residuals of the Mbh - o" 
correlation against Rf. that show a weak, but significant, correla- 
tion (reproduced here in Fig.ll6>. The new measurement of the 
black hole mass in NGC 3998 supports this idea. In fact, NGC 
3998 has one of the smallest values of R^ among galaxies with 
measured Mbh (0.85 kpc) and it sh ows a negative res idual from 
the Mbh - cr correlation. Recentlv iCaoetti et al.l ( 120051) found a 
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Fig. 11. Best fit values of Ty and Mbh obtained by varying the 
disk inclination / (upper and middle panels). The corresponding 
are shown in the bottom panel. The lines in the upper and mid- 
dle panels reproduce the oc 1/sin^ / dependence of Mbh and Ty. 
Inclinations larger than 70° are excluded at a confidence level 
of 2cr (shaded region in the right side of the diagrams). Models 
with inclinations smaller than 27° are excluded since they corre- 
spond to unacceptably high values of Yy (shaded region in the 
left side of the diagrams, bounded by Ty = 7.87 as expected for 
a Salpeter IMF). 
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Fig. 13. Line width distribution expected from the overall best 
fit model (solid line) at / - 30°, and derived (dashed line) by 
including the observed line widths in the fitting procedure. From 
upper to bottom panel; N2, Nl, NUC, SI and S2. 
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Fig. 14. Mbh vs bulge mass from lMarconi & Hun3 (l2003h with 
the best fit obtained from a bisector linear regression analysis 
(solid line). The position of NGC 3998 (empty circle) is indi- 
cated by an arrow. 



Fig. 12. Xi contours at varying Ty and Mbh at an inclination of 
/ = 30°. Contours are plotted for confidence levels of 1, 2 and 3 
cr. The plus sign marks the overall best fit. 



similar result, but in the opposed sense, considering the Seyfert 
galaxy NGC 5252: a large effective radius (9.7 kpc) corresponds 
in this galaxy to a large positive residual. This confirms that a 
combination of both cr and 7?e is necessary to drive the correla- 
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Log a [km s ^] 

Fig. 15. Mbh vs ce ntral stellar velocity dispersion from 
iTremaine et ajj J2002h (upper panel) and from lFerrarese & FordI 
(2005) (lower panel). Solid lines indicate the best fits for the two 
correlations. The position of NGC 3998 is shown by an empty 
circle. 

tions between Mbh and other bulge properties, an indication for 
the presence of a black holes "fundamental plane". 

6. Summary and conclusions 

We have presented results from a gas kinematics study in the nu- 
cleus of the nearby SO active galaxy NGC 3998. The analyzed 
data were retrieved from archival HST/STIS long-slit spectra. 
We performed an analysis of the Ha, [N 11]/1/16548,6583 and 
[SII]/l/l6716,6731 emission lines profiles to derive the map of 
the gas velocity field. The nuclear velocity curves show a gen- 
eral reflection symmetry and are consistent with the presence 
of gas in regular rotation. We used our modeling code to fit 
the observed Ha surface brightness distribution and velocity 
curve. The dynamics of the rotating gas can be accurately repro- 
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Fig. 16. Residuals from the Mbh vs cr correlation (in the 
Tremaine's et al. form) reported against the galaxy's eff'ective 
radius R^. The position of NGC 3998 is indicated by an empty 
circle. 



duced by motions in a thin disk when a compact dark mass of 
Mbh = 2.7;^2o ^ 10'^ M©, very likely a supermassive black hole, 
is added to the stellar mass component. This result is also sup- 
ported by the satisfactory good match of observed and model 
line width distributions. Furthermore, the black hole in NGC 
3998 has a sphere of influence radius of ~ 13 pc (0'.'16). At the 
high HST spatial resolution, this value of R^ph implies a resolved 
BH sphere of influence for our Mbh determination (27?sph/^res - 
3.2). 

For what concerns the connections of this BH mass estimate 
with the properties of the host galaxy, the Mbh value for NGC 
3998 is in excellent agreement (within a factor of 2) with the 
Mbh - Mbui correlation between BH and host bulge mass. The 
black hole mass predicted by the Mbh - crbui correlation is a fac- 
tor of 2.7 larger than our measure, adopting the relation found 
by Tremaine et al. ( 200 ^, or a factor o f 7 usin g the most recent 
parameterization by F errarese & FordI (l2005l) . However, NGC 
3998 cannot be considered as an outlier from the Mbh vs. cr cor- 
relation considering the errors in the estimates of Mbh and osmr, 
and both the uncertainties in the determination of the Mbh - o" 
relation as well as its scatter. 

Nonetheless, the lower-than-expected value for the mass of 
the black hole hosted by NGC 3998 strengthens the presence of a 
connection between the residuals from the Mbh - o" relation and 
the galaxy's effective radius. In fact, NGC 3998 has one of the 
smallest values of R^ among galaxies with measured Mbh and it 
shows a negative residual. We also recently showed that the op- 
posite is true for the Seyfert galaxy NGC 5252: a large eff'ective 
radius corresponds in this galaxy to a large positive residual. 

Apparently only with a combination of both cr and 7?e it 
is possible to account for the correlations between Mbh and 
other bulge properties, an indication for the presence of a black 
holes "fundamental plane". Clearly, only by further increasing 
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the number of direct black hole measurements it will be possible 
to base these conclusions on a stronger statistical foundation. 
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Appendix A: Fitting parameters for STIS spectra. 
List of Objects 

INGC 3998f on page 1 
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Table A.l. Fitting parameters for the narrow emission line spectra at each location (x coord.) along the slit. The relative flux of the 
[Nil] lines was held fixed to 0.334, so only the [NII]/16583 is hsted. From upper to bottom table: NUC, Nl, SI, N2 and S2 positions. 



X 


V 


dV 


FWHM 


dFWHM 


I-iHa) 


dE{Ha) 




dUiNII) 


msii) 


dLiSII) 




dI.{SII) 
















(/16583) 


(i6583) 


(i6716) 


(^6716) 


(^6731) 


(A6731) 


591-595 


1204 


23 


590 


35 


61 


7 


77 


6 


25 


5 


24 


6 


596 


1163 


21 


561 


39 


135 


16 


161 


13 


42 


10 


63 


11 


597 


1163 


16 


495 


38 


143 


24 


198 


24 


89 


9 


77 


10 


598 


1139 


11 


541 


29 


208 


26 


288 


25 


142 


10 


150 


11 


599 


1157 


14 


750 


37 


934 


71 


1060 


72 


262 


25 


347 


30 


600 


1178 


10 


814 


30 


1853 


109 


2257 


113 


577 


32 


701 


40 


601 


1117 


13 


1044 


28 


6873 


230 


6210 


168 


883 


78 


1647 


93 


602 


867 


8 


1085 


23 


5315 


107 


4187 


94 


766 


31 


1096 


37 


603 


837 


8 


646 


21 


1081 


55 


1140 


51 


311 


16 


387 


20 


604 


882 


6 


422 


16 


269 


25 


325 


20 


158 


8 


156 


10 


605 


905 


7 


391 


16 


169 


18 


224 


15 


117 


7 


129 


9 


606 


917 


8 


434 


18 


200 


14 


220 


12 


107 


8 


110 


9 


607 


917 


10 


431 


21 


134 


11 


150 


10 


79 


6 


80 


8 


608 


894 


14 


407 


32 


82 


9 


93 


6 


44 


7 


41 


8 


609 


870 


11 


331 


23 


67 


8 


78 


8 


39 


5 


37 


6 


610 


842 


14 


333 


32 


51 


8 


54 


7 


39 


5 


16 


6 


592-596 


1172 


15 


548 


29 


83 


6 


387 


5 


167 


4 


152 


5 


597 


1111 


24 


598 


46 


133 


17 


161 


13 


43 


11 


59 


12 


598 


1118 


12 


594 


24 


257 


15 


268 


12 


65 


9 


92 


10 


599 


1159 


10 


569 


22 


356 


22 


421 


22 


131 


9 


184 


11 


600 


1163 


9 


516 


24 


532 


51 


553 


41 


172 


15 


216 


18 


601 


1168 


8 


521 


21 


551 


49 


622 


40 


215 


13 


251 


15 


602 


1119 


12 


615 


32 


576 


70 


585 


53 


214 


12 


211 


14 


603 


1018 


20 


592 


47 


292 


80 


396 


57 


195 


12 


194 


15 


604 


946 


7 


428 


16 


197 


18 


251 


16 


124 


7 


145 


9 


605 


948 


6 


351 


15 


154 


14 


162 


12 


99 


6 


86 


7 


606 


942 


7 


359 


18 


137 


11 


162 


10 


77 


7 


77 


8 


607 


954 


7 


304 


18 


88 


8 


99 


7 


62 


5 


64 


6 


608 


935 


9 


289 


27 


57 


8 


80 


8 


35 


5 


48 


6 


609 


923 


14 


325 


28 


57 


8 


61 


7 


29 


5 


10 


5 


610 


929 


16 


344 


31 


36 


6 


46 


6 


21 


4 


28 


5 


592-596 


1130 


11 


453 


27 


72 


6 


76 


5 


39 


4 


32 


4 


597 


1094 


12 


450 


32 


145 


13 


140 


20 


54 


8 


65 


9 


598 


1075 


6 


348 


17 


154 


11 


183 


9 


56 


7 


88 


8 


599 


1062 


9 


409 


23 


161 


24 


231 


20 


98 


8 


125 


10 


600 


1050 


10 


622 


27 


469 


29 


476 


28 


134 


11 


192 


13 


601 


1001 


11 


730 


35 


580 


32 


657 


34 


149 


12 


222 


15 


602 


914 


11 


718 


33 


633 


35 


629 


34 


144 


13 


226 


16 


603 


888 


9 


481 


24 


373 


36 


328 


28 


123 


9 


116 


10 


604 


885 


8 


406 


21 


196 


19 


214 


17 


108 


7 


77 


8 


605 


879 


10 


431 


25 


160 


17 


170 


16 


87 


8 


62 


9 


606-610 


867 


8 


458 


17 


82 


5 


92 


5 


45 


3 


46 


4 


592-595 


1099 


22 


515 


45 


38 


6 


58 


5 


12 


4 


15 


5 


596-598 


1101 


17 


512 


37 


86 


9 


94 


7 


34 


6 


45 


6 


599-600 


1090 


14 


489 


35 


106 


11 


134 


9 


33 


6 


60 


8 


601-602 


1077 


11 


498 


24 


123 


10 


166 


9 


48 


6 


63 


7 


605-606 


994 


9 


360 


23 


78 


8 


110 


7 


38 


5 


48 


6 


607-610 


980 


11 


372 


26 


46 


5 


53 


5 


23 


3 


31 


4 


592-594 


1095 


11 


342 


27 


49 


6 


50 


4 


24 


3 


28 


4 


595-596 


1100 


15 


468 


30 


81 


9 


85 


7 


38 


5 


39 


6 


597-598 


1079 


18 


519 


41 


85 


10 


95 


8 


35 


6 


24 


7 


599-600 


1018 


19 


549 


37 


96 


11 


115 


9 


38 


6 


41 


7 


601-602 


1017 


20 


541 


37 


90 


10 


97 


8 


30 


5 


37 


7 


603-605 


921 


21 


548 


41 


86 


9 


81 


7 


29 


5 


27 


6 


606-610 


811 


15 


383 


33 


34 


5 


42 


4 


16 


3 


18 


3 



